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ABSTRACT 

From the parallaxes and proper motions of a kinematically unbiased subsample of the 
Hipparcos catalogue we have redetermined as a function of colour the kinematics of 
main-sequence stars. 

Whereas the radial and vertical components of the mean heliocentric velocity of 
stars show no trend with colour, the component in the direction of galactic rotation 
nicely follows the asymmetric drift relation, except for stars bluer than B—V = 0.1 mag. 
Extrapolating to zero dispersion yields for the velocity of the Sun w.r.t. the local 
standard of rest (LSR) in kms" 1 : U = 10.00±0.36 (radially inwards), V Q = 5.23±0.62 
(in direction of galactic rotation), and Wo = 7.17 ± 0.38 (vertically upwards). 

Parenago's discontinuity is beautifully visible in a plot of velocity dispersion 
against colour: the dispersion, which is essentially constant for late spectral types, 
decreases towards early spectral types blueward of B — V ~ 0.62 mag. 

We determine the velocity-dispersion tensor <r 2 as function of colour. The mixed 
moments involving vertical motion are zero within the errors, while crjL i s non-zero at 
about (lOkms -1 ) 2 independent of colour. The resulting vertex deviations are about 
20 degrees for early-type stars and 10 ± 4 degrees for old-disc stars. The persistence 
of the vertex deviation to late-type stars implies that the Galactic potential is sig- 
nificantly non-axisymmetric at the solar radius. If spiral arms are responsible for the 
non-axisymmetry, they cannot be tightly wound. 

Except for stars bluer than i3-V = 0.1mag the ratios of the principal velocity 
dispersions are given by <r\ : a 2 ■ 03 ~ 2.2:1.4:1, while the absolute values increase with 
colour from <j\ rs 20 kms -1 at B-V = 0.2 mag to <j\ rs 38 kms -1 at Parenago's disconti- 
nuity and beyond. These ratios imply significant heating of the disc by spiral structure 
and that R^/Rd — 3 to 3.5, where Rd is the scale length of the disc. 

Key words: Stars: kinematics - Galaxy: fundamental parameters - Galaxy: kine- 
matics and dynamics - Galaxy: solar neighborhood - Galaxy: structure 



1 INTRODUCTION 



The kinematics of stars near the Sun has long been known to 
provide crucial information regarding both the structure and 
the evolution of the Milky Way. Karl Schwarzschild (1908) 
already interpreted the distribution of random velocities as 
forming a triaxial 'velocity ellipsoid', which Oort, B. Lind- 
blad and Stromberg were able to relate to the large-scale 
structure of the disc. In the early 1950s Parenago (1950), 
Nancy Roman (1950, 1952) and others pointed out that 
stellar kinematics varies systematically with stellar type, in 
the sense that groups of stars that are on average younger 
have smaller velocity dispersions and larger mean Galac- 
tic rotation velocities than older stellar groups. Spitzer & 
Schwarzschild (1953), Barbanis & Woltjer (1967), and Wie- 
len (1977) explained these correlations in terms of the diffu- 
sion of stars through phase space as the Galactic disc ages 



- for recent studies of these processes see Binney & Lacey 
(1988) and Jenkins (1992). 

The Hipparcos Catalogue (ESA 1997) provides an im- 



portant opportunity to re-examine the fundamental data of 
solar-neighbourhood kinematics by providing the first all- 
sky catalogue of absolute parallaxes and proper motions. 
From the Hipparcos Catalogue we can, moreover, extract 
samples that are completely free of the kinematic biases 
that have plagued similar studies in the past. Binney et al. 
(1997) obtained some preliminary results from a sample of 
5610 stars around the south celestial pole. Here we extend 
these results in three ways: 

1. We have joined the relatively deep sample of Binney et 
al. (1997) to a sample that covers the whole sky but has 
a significantly brighter limiting magnitude. The combined 
sample gives significantly better statistics for earlier spectral 
types. 
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Figure 1. Hcrtzsprung-Russcll diagram (Mh p is the absolute 
magnitude in Hipparcos's own passband) of the 18 860 single Hip- 
parcos stars with relative parallax errors less than 10 per cent. 
The lines are used to select the main sequence and have 16 054 
stars between them. 



2. Since we now have some all-sky coverage we are able to 
impose a stricter limit on parallax errors: we use only stars 
with ov/tt < 0.1, whereas Binney et al. accepted stars with 
(j^/tt < 0.2. 

3. Whereas Binney et al. only reported the solar motion 
and a single estimate of random velocity for each spectral 
type, we derive the lengths of the axes of the velocity ellip- 
soid and analyze its orientation. 

Section 2 describes the new sample. Section 3 explains 
how we have analyzed the sample. Section 4 gives the results 
and Section 5 discusses them. 



2 THE SAMPLE 

Our sample has to meet several criteria if it is accurately 
to determine local stellar kinematics as a function of stellar 
age. First, it must be kinematically unbiased, i.e., faithfully 
represent at each point in the Hertzsprung-Russell diagram 
the kinematics of all nearby stars. Second, it should be con- 
fined to main-sequence stars, since only for these is there a 
one-to-one relation between age and either colour or absolute 
lumino sity. Third, it should be based on accurate astrome- 



try, anc , fourth, it should contain no multiple stars because 




Figure 2. Hertzsprung-Russell diagram for the stars in our kine- 
matical unbiased sample of 11 865 single main-sequence stars. 



the Hipparcos catalog. The Hertzsprung-Russell diagram of 
the selected stars is shown in Figure [j], which also shows the 
division lines that we used to select main-sequence stars; 
there are 16 054 stars between these lines. However, given 
the very heterogeneous nature of the Hipparcos catalog, they 
almost certainly do not comprise a kinematically unbiased 
sample. The only reliable way to extract such a sample is to 
isolate a magnitude-limited subsample. 

According to its Volume I, p. 4, the Hipparcos Cata- 
logue comprises some 60 000 objects complete to about 7.3 
to 9 mag depending on galactic latitude, b, and stellar type. 
Actually, however, 59 stars brighter than Vt = 7.0 (the sub- 
script T denotes the passbands used by Tycho, the starmap- 
per of Hipparcos) are absent from the Catalogue (see the 
t able on p. 142 of Volume I) . We used the Tycho Catalogue 
(ESA 1997), which is essentially complete to 11 mag, to con- 
struct an almost complete subsample of Hipparcos stars. For 
each of 16 x 16 x 10 uniformly spaced bins in sin b, galac- 
tic longitude, £, and Bt — Vt, we took all Hipparcos stars 
brighter (in V) than the second brightest star per bin that 
is contained in the Tycho but not the Hipparcos Catalogue. 
This gives a subsample with 95 per cent completeness and 
47 558 stars (compared to 60 000 claimed), of which 10 706 
have o^/tt < 0.1 and lie between the lines on Figure [j]. 

Another kinematically unbiased subsample is given by 
Hipparcos Proposal 018: all st ars spectrally classified in the 
Mich i gan Catalo gue by 1982 ([Houk fc Cowley 1975|; |Houk 



their kinematics contains additional motions. 

We ensured satisfaction of the last two criteria by taking 
only single stars with relative parallax errors smaller than 
10 per cent. This selection criterion automatically excludes 
remote stars, and leaves just 18 860 of the 118 218 stars in 



1975; Houk 1982) that, judged from their spectral classifica- 
tion, should be within 80 pc from the Sun - this 'SCP' sam- 
ple is the one analyzed by Binney et al. (1997). It contains 
6845 stars south of S = —28 degree (covering 26.5 per cent 
of the sky), of which 3172 are single and have a^/n < 0.1 
and lie between the lines on Figure oj. 
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Whereas most of the stars in the all-sky sample are 
of early type and intrinsically luminous, the SCP sample 
contains many late-type dwarfs. The union of the samples 
is kinematically unbiased and has 11865 entries. It is this 
combined sample that we analyze below. Figure |^ shows its 
Hertzsprung-Russell diagram. 



3 ANALYSIS TECHNIQUE 

Hipparcos provides us with parallax measurements of un- 
precedented accuracy, and proper motions whose accuracy 
is only comparable to that of ground-based astrometry on 
account of short duration (3.3 yr) of the Hipparcos mission. 
Unfortunately, the Hipparcos astrometry mission was not 
complemented by a programme to measure the radial velo- 
cities of the same stars, so we do not know the space ve- 
locities of most sample stars. As Binney et al. (1997) have 
demonstrated, the stars with known radial velocities form a 
kinematically biased subsample. Hence we have to discard 
the known radial velocities and work with the Hipparcos 
data alone. 



3.1 The projection equation 

The velocity of a star relative to the Sun can be divided 
into three components: (i) the peculiar velocity of the star, 
(ii) that of the Sun, and (iii) the contribution from Galactic 
rotation. Here we are dealing with the first two parts, and 
correct for Galactic rotation before further analysis. Given 
observed value fi e ° and fi)° bs ', the corrected values are 



/i f : 
fib : 



(obs) 



Acos(2l) - B 



: M 6 



(obs) + Asin(2^) cos b sin b. 



(1) 



We used the values of A and B derived by Feast & Whitelock 
(1997) from Hipparcos Cepheids: A = 14.82 km s -1 kpc -1 
and B = -12.37 km s _1 kpc -1 . 

Let us introduce a Cartesian coordinate system such 
that x points towards the Galactic centre, y in the direction 
of Galactic rotation, and z towards the north Galactic pole. 
We consider a star with parallax 7r and heliocentric velocity 
v, and define the star's proper-motion velocity p to be 



P : - 



- sin £ cos b — cos I sin b fit, 
cos £ cos b fie — sin £ sin b fit 
cos b fib 



Then p and v are related by 
p = A ■ v, 

where the matrix A is defined by 
A = I — r <g> f, 



(2) 



(3) 



(4) 



with r the unit vector to the star. A projects velocities onto 
the celestial sphere. It is symmetric and, as every projection 
operator, obeys A 2 = A and is singular. Hence we cannot 
invert (B|); we needed the radial velocity to recover v. 



3.2 The mean motion 

The solar motion relative to any given group of stars is sim- 
ply minus the mean motion of that group with respect to 



the Sun, vq = — (v) . For a kinematically unbiased sample 
of nearby stars we can safely assume that the positions on 
the sky r are uncorrelated with the velocities v. With this 
assumption, taking the sample mean of Equation yields 
(p) — (A) ■ (v), which can be inverted to give 

- v e = (v) = (A)' 1 ■ (p) . (5) 

It is easy to show that for a group of stars that is isotropically 
distributed over a celestial hemisphere (A) = 2 /:iI, where 
the factor 2 /s arises from the fact that we only know two of 
the three velocity components for each star. Let us denote 
the motions relative to the mean by 



v' = v — (v) 
p' = p — A • (v) 

and consider the quantity 

s 2 = (\p'\ 2 )- 



(6) 



(7) 



It can be shown that the choice of (v) given by (||) minimizes 
S 2 , and (with Equations ^| and ^) that 



s a = (M a >-((«'.f) a >. 



(8) 



Hence S 2 is a measure for the velocity dispersion of the 
group. 

The variance of (v) can be estimated to be 

->2 



V ({v}) = N' 1 {A}' 1 S 2 



(9) 



where N is the number of stars in the sample, while an 
estimate for the variance of S 2 is 



v(s 2 ) = n- 1 (OpTMIp'I 2 ) 2 ) 



(10) 



These error estimates only account for the Poisson noise due 
to the finite number of stars, and measurement uncertainties 
in p have been neglected. 

3.3 The velocity dispersion tensor 

Similar to the mean velocity, the second-order moments of 
the space velocities can be inferred. From Equation Q we 
can obtain 

p ®p = (A - v) ® (A ■ v') 

= (A <g> A) • (v' <g> v). (11) 

The velocity dispersion tensor cr 2 is the sample average of 
(i/cg>i/), so after averaging (11) we can solve for cr 2 , exactly 
as we did for (v) above, and the same holds, in principle, for 
the velocity moments of higher orders. 

However, if we were to solve the expectation of Equa- 
tion (11) for cr 2 , we would in general find that the recovered 
tensor cr 2 was not symmetric, because (AjjAki) (AkjAu). 
Since we know a priori that cr 2 is a symmetric tensor, with 
six rather than nine independent elements, it is advanta- 
geous to impose this symmetry by solving 



(PiPk) = I X] (AijAki + AkjAu) crji. 



(12) 



For the computer this is more conveniently written as (it) = 
B ■ (s) , where (u) and (s) are vectors containing the six inde- 
pendent components of (p'ip'j) and cry, respectively, whereas 
B is the corresponding 6x6 matrix. With this notation the 
dispersion tensor is estimated via 
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Figure 3. The components U, V, and W of the solar motion w.r.t. stars with different colour B-V. Also shown is the variation of the 
dispersion S with colour. 



(s) = (B)- 1 ■ («> , (13) 

while an estimate for its variance is given by 

y« S )) = 7V- 1 (B)- 1 (| W -B.( S )| 2 ). (14) 

Somewhat surprisingly, for an isotropic sample the matrix 
(B) _1 is not diagonal, in contrast to (A) -1 . It turns out 
that the diagonal terms of cr 2 are coupled. Quantitatively, 



x, y, and z directions as defined in Section 3.1. U and W 




(PxPx) 
(PyPy) 

ip'zPz), 



(15) 



while all other off-diagonal terms of (B) 1 average to zero 
for an isotropic sample. 



4 RESULTS 

We will now apply the formulae from the last section to our 
sample. The analysis technique relies heavily on taking aver- 
ages of the various observed kinematic quantities. Unfortu- 
nately, averages are sensitive to outliers such as halo stars, 
and rejection of them is an important issue. Here, we will 
use an iterative method: in the solution for the various quan- 
tities we rejected stars that contribute to S 2 (Equation ^) 
by more than k times the value of S 2 obtained in the pre- 
vious iteration. Convergence is usually obtained after a few 
iterations (or never if k is too small). The results prove to 
be insensitive to k in the range of 3 to 4; the results below 
were obtained with = 4. 



4.1 The solar motion and Parenago's discontinuity 

We have binned the stars in B-V , such that the bins are no 
smaller than 0.02 mag and contain no less than 500 stars. 
For each bin, Figure [3| plots the solar motion with respect 
to the stars in the bin, vq = — (v), and S, which is a mea- 
sure for the bin's velocity dispersion, versus the mean colour. 
As usual, U, V, and W denote the components of vq in the 



do not vary significantly between bins, while both V and 
S increase systematically from early to late spectral types. 
The points in S display very beautifully Parenago's (1950) 
discontinuity: around B-V ~ 0.62 mag there is an abrupt 
change in gradient from a strongly positive value to about 
zero. The same discontinuity is visible, though less clearly, 
in the data for V. Parenago's discontinuity is thought to 
arise from the fact that the mean age of stars decreases as 
one moves blueward of the discontinuity, while it is inde- 
pendent of colour redward of the discontinuity: scattering 
processes cause the random velocities of stars to increase 
steadily with age (e.g. Jenkins 1992). Hence velocity disper- 
sion reflects age, and decreases as one moves blueward from 
the discontinuity through ever younger stellar groups, while 
remaining constant with mean age redward of the disconti- 
nuity. The discontinuity itself should occur at the colour for 
which the main-sequence lifetime of a star equals the age of 
the Galactic disc. However, since stars change colour during 
their life on the main sequence, detailed modelling of stellar 
populations is necessary, to infer the age of the stellar disc 
from this datum (Binney & Dehnen, in preparation). 



4.2 The Sun's velocity w.r.t. the LSR 

Figure m is a plot of U, V and W versus S 2 . For S > 
15kms , this clearly shows the linear dependence of V 
on S 2 that is predicted by Stromberg's asymmetric drift 
equation [e.g. Equation (4-34) of Binney & Tremaine, 1987 
(hereafter BT)]. That is, V increases systematically with S 2 
because the larger a stellar group's velocity dispersion is, 
the more slowly it rotates about the Galactic centre and the 
faster the Sun moves with respect to its lagging frame. 

For very early- type stars with B-V SS 0.1 mag and/ or 
S is 15kms _1 , the V-component of (v) decreases with in- 
creasing S, colour, and hence age contradicting the expla- 
nation given in the last paragraphs. However, the stars con- 
cerned are very young, and there are several possibilities for 
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Figure 4. The dependence of U, V, and W on S 2 . The dotted lines correspond to the linear relation fitted (V) or the mean values (U 
and W) for stars bluer than B-V=0. 



them not to follow the general trend. First, because of their 
youth these stars are unlikely to constitute a kinematically 
well-mixed sample, rather they move close to the orbit of 
their parent cloud; many will belong to a handful of moving 
groups. Second, Stromberg's asymmetric drift relation pre- 
dicts a linear relation between V and S 2 only if both the 
shape of the velocity ellipsoid, i.e. the ratios of the eigen- 
values of er, and the radial density gradient are independent 
of S. Young stars probably violate these assumptions, espe- 
cially the latter one. 

The velocity (Uo, Vq, Wo) of the Sun with respect to the 
local standard of rest (the velocity of the closed orbit in the 
plane that passes through the location of the Sun) may be 
read off from Figure ^ by extrapolating any trends of U, V 
and W with S 2 back to S = 0. Ignoring stars blueward of 
B-V = mag we find 



Uo = 10.00 ± 0.36 kms 
Vb = 5.23 ± 0.62 kms 
W = 7.17 ±0.38 kms 



(16) 



When we use this value of the solar motion to calculate 
values of (v y ) relative to the LSR for our stellar groups, 
Stromberg's asymmetric drift equation is found to be 



<«y) = 



(17) 



where k — 80 ± 5 kms 1 and we have used S 2 ~ l.lAo 2 . 



4.3 The velocity dispersion tensor 

Once the Sun's velocity vq with respect to the LSR is 
known, we can determine proper motions p' with respect to 
the LSR by replacing (v) by —v@ in Equation (^). Then for 
any stellar group we can determine a 2 from these values of 
p' as described in Section 3.3. We divided the 11 865 stars of 



our sample in nine bins in B-V with equal numbers of stars 
in each bin. The results are displayed in Figure ^| and Ta- 
ble 0. The errors correspond to the 15.7 and 84.3 percentiles 
(la error) and have been evaluated assuming a multivariate 



Table 1. Eigenvalues of a 2 for the nine colour bins and all stars 
beyond Parenago's discontinuity (last row). 



bin (B-V) 



mi ii, max 



Cl 



(Tl/<T2 0"l/f"3 



-0.238 0.139 

0.139 0.309 

0.309 0.412 

0.412 0.472 

0.472 0.525 

0.525 0.582 

0.582 0.641 

0.641 0.719 

0.719 1.543 



14 40+ ' 49 
20.23+° 43 
22.40l -j« 
26.331°;*° 
30 45+ - 96 

1.08 
0.75 
+ 1.37 
0.94 
1.19 
0.85 
1.40 
0.93 



1.52 
2.10 



+0.16 

0.14 
+0.13 
0.28 
+ 0.13 



33.02 
37.73 
38.23 
37.28 



L88 -0.20 

1 65+ 012 
1-OO-0.15 

1 66+ 013 

± ' ol -0.12 

1 60+ 07 

1 ' ou -0.18 
1 59+ - 08 
1 43+0.12 



2 62+ 91 
z -° -0.28 

2 50+ ' 81 
^• ou _0.10 

9 qn+0-65 

^• oa -0.14 

+0.60 

-0.14 

+0.76 

0.18 

+0.63 

0.19 

+0.47 

0.04 

+0.38 

0.06 

2 04+ ' 60 



22.; 



2.15 
2.27 
2.18 
1.77 
1.83 



3o.3t* : ; 

+2.8 
-3.0 

19.813 4 

10.21^? 
in o+ 5 - 6 

'• D -5.5 
13 1+ 67 



0.620 1.543 37.97+ 



0.81 



1.52 



+0.08 
0.09 



3.9 
4.1 



fli C2> ""3 are the roots of the largest, middle, and smallest eigen- 
value of the velocity dispersion tensor a 2 . l v is the vertex devia- 
tion (Equation |ri|). Units are mag, kms -1 , and degree for B-V, 
<Ji, and l v respectively. The errors given correspond to the 15.7 
and 84.3 percentiles, i.e. ltr error. 

Gaussian distribution in the a 2 j with variance evaluated via 
Equation (|l4|). In the upper panel of Figure ^|the three di- 
agonal velocity dispersions and (v^) are plotted versus B-V. 
Parenago's discontinuity is visible in all three era and (v$), 
though, due to the larger bin size, less clearly than in Fig- 
ure H above. The ordering between the diagonal components 
of <j is the same for all colour bins: a xx > a yy > a zz . 
The lower three panels of Figure |^ show 



ff« = sign( 



(18) 



for the mixed components of the velocity dispersion tensor 
cr 2 . Evidently, the mixed moments involving vertical mo- 
tions vanish within their errors. This is to be expected for 
essentially all possible dynamical configurations of the Milky 
Way. On the other hand, the mixed dispersion in the plane, 
a 2 y differs significantly from zero, which is not allowed in 
a well-mixed axisymmetric Milky Way. Thus, the principal 
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Figure 5. Velocity dispersions for stars in different colour bins. 
The top panel shows the mean rotation velocity (negative values 
imply lagging w.r.t. LSR) and the three main velocity dispersions. 
In the three bottom panels a+ = sign(cr? j . ) | c?^ | 1 / 2 is plotted for 
the mixed components of the tensor a 2 .. 



axes of cr 2 are not aligned with our Cartesian coordinate 
frame. We diagonalized the tensor cr 2 to obtain its eigenval- 
ues erf. The square root of the largest of these as well as the 
ratios to the smaller ones are given in Table |l|. The ratio 
eri/<72 — 1.6, whereas cti/ctz ~ 2.2 with a trend for smaller 
values at redder colours. 

Also given in Table [j] is the "vertex deviation", com- 
monly used to parameterize the deviation from dynamical 
symmetry. This is defined to be 



2at 



(19) 



and is the angle by which one has to rotate our Cartesian 
coordinate system around its z axis such that the result- 
ing velocity dispersion tensor is diagonal in the v x v y plane. 
Hence £ v is the Galactic longitude of the direction of a 2 , the 
largest eigenvalue of the velocity dispersion tensor. Fig. [] 
shows £ v as a function of colour. Blueward of Parenago's 
discontinuity there is a clear trend of £ v decreasing with B- 
V. Redward of Parenago's discontinuity stellar kinematics is 



•10 



, 30 



20 



:o 



0.2 0.4 0.6 

B-V [mag] 



0.8 1 



Figure 6. The vertex deviation £ v versus B-V colour. The error 
bars correspond to the 15.7 and 84.3 percentile (i.e. lc error) and 
have been obtained assuming a multivariate Gaussian distribution 
in the u 2 ^ with variance evaluated via Equation (p^). 



independ of B-V and it makes sense to group all 3093 stars 
into a single bin. The results are given in the last row of 
Table [j]. The hypothesis that £ v < for this gsoup of stars 
is excluded at the 99 per cent confidence level.El 



5 DISUSSION 

We have analyzed a kinematically unbiased sample of nearly 
12 000 main-sequence stars for which the Hipparcos satellite 
measured parallaxes that have errors smaller than 10 per 
cent. From this sample we have redetermined most of the 
key kinematic parameters of the solar neighbourhood with 
unprecedented accuracy. 

On account of the unusually small accidental and sys- 
tematic errors of the Hipparcos data, our errors are domi- 
nated by Poisson noise; we are characterizing the distribu- 
tion of stars in three-dimensional velocity space from sam- 
ples of typically 500 stars. From this fact it follows that our 
error estimates may be considered reliable. 

Equation ( |l6|) gives the velocity of the Sun with re- 
s pect to the LSR. In common with several recent s tudies 
ilavor 1974 IComez fc Messenier 1977[ |Oblak 1983 



Bien 



aymc &r Sechaud 1997| ) this has a smaller azimuthal compo- 



nent (Vo = 5.2 ± 0.6 km s -1 ) than Delhaye's classical value, 
12kms . A contributory factor to this discrepancy is prob- 
ably the fact that stars bluer than B-V — have anomalous 
values of both (v x ) and {v y }, presumably because they are 
not yet well mixed. We have ignored such stars in our de- 
termination of the Sun's motion with respect to the LSR. 

Equation (^) quantifies the asymmetric drift. As is well 
known, the coefficient k that occurs in it constrains the ra- 
dial gradients of stellar density and the velocity-dispersion 
tensor, cr 2 . Quantitatively, if i?a is the scale-length of the 
disc, and v the mid-plane stellar density, with Equation (4- 
34) of BT we may infer from Equation ([17|) that 



Ro 



1 d In ual x 

2 dhxR 



fa+(i-£p 



(20) 



: 3.0 + 0.4. 



It is important to bear in mind that the errors in £ v are not 
distributed Gaussianly 
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Here the first equality involves the assumption that a 2 x oc v, 
which is based on the observation that disc scale heights 
vary little with radius, and the term that follows the symbol 
+ may or may not be required, depending on whether the 
longest axis of cr 2 points to the galactic centre or remains 
horizontal as one moves above the plane [cf. §4.2.1 (a) of BT]. 
The numerical values given are based on Table 1, Ro = 8 kpc 
and the Feast & Whitelock values of the Oort constants. 
The relatively short scale-length implied by Equation (20) 
i s consistent with recent inrared-based studies of t he Galaxy 



( |Kent, Dame fc Fazio 199jJ ; |Malhotra et al. 1997|) 



The shape of the velocity ellipsoid does not vary sig- 
nificantly with B-V; to within the errors its axis ratios 
are constant at 01/03 ~ 2.2 and <7i/<72 — 1.6. The ratio 
<jy y /(Jxx — 0.4 is intimately connected with the values of 
the Oort constants, being naively given by Oort's equation 
—B /(A — B) ~ 0.45. Unfortunately, Oort' s relation is s ignif- 
icantly in error for typical star samples (Binney 1986) but 
Cuddeford & Binney (1994) argue that —B/(A — B) can be 
accurately determined from <Jy y /a xx if modified rather than 
normal moments are employed. We hope soon to evaluate 
these modified moments from the present sample. 

The ratio 0-1/0-3 constrains the nature of the scattering 
process that is responsible for the increase in stellar velocity 
dispersion with age (e.g. Jenkins 1992). The value deter- 
mined here is consistent with previous values, and implies 
that both spiral structure and scattering by molecular clouds 
contribute significantly to heating of the local disc. We find 
a hint of smaller values of 01/03 for redder colours, which 
is the trend expected if heating by spiral structure is less 
efficient relative to cloud-scattering for dynamically hotter 
stellar populations. 

We have shown that one principal axis of the velocity 
ellipsoid coincides with the direction b = 90deg. The angle 
by which the longest axis of the velocity ellipsoid deviates 
from the direction to the Galactic centre, the vertex devia- 
tion £ v , decreases from large values (£ v > 20deg) for early 
types up to the colour of Parenago's discontnuity and is then 
consistent with being constant at £ v ~ lOdeg. 

Two factors probably contribute to the vertex devia- 
tion: (i) the likelihood that a significant fraction of young 
stars belong to a small number of moving groups, and (ii) 
any large-scale non-axisymmetric component in the Galac- 
tic potential, such as would be generated by either a bar 
or spiral structure. Spiral arms will contribute to the vertex 
deviations of all stellar groups, but their contribution will 
be largest for the populations with the smallest velocity dis- 
persions because stars with epicycle amplitudes comparable 
to or larger than the interarm spacing will respond weakly 
to the spiral's potential. Hence the more tightly wound the 
arms are, the more strongly the effect will be confined to 
the earlier spectral types. The fact that we see a signifi- 
cant vertex deviation to the latest spectral types, implies 
that the causative agent is either a bar or rather open spiral 
arms. To quantify this point, let X be the epicycle amplitude 
of a typical late-type star. Then the star's radius is R(t) = 
R g + Xcos(Kt), where k = 2V B 2 - AB ~ 36.7 km s" 1 kpc" 1 
is the epicycle frequency. To a reasonable approximation we 



may equate the time average of R' 



to u\ 



(38 kms -1 ) 2 



the orbit of such a star. It is clearly important to quantify 
this interesting conclusion more precisely. 

Finally it seems worth remarking that in characteriz- 
ing the distribution of stars in velocity space by its first few 
moments we in no way imply that the distribution is similar 
to Schwarzschild's ellipsoidal distribution, as given, for ex- 
ample, by equation (7-91) of BT; the velocity ellipsoid is a 
formal construct which need have no physical counterpart. 
From these data it is possible to map the distribution of 
stars in velocity space and discover how closely it resembles 
Schwarzschild's paradigm. But that topic is reserved for a 
future paper (Dehnen, in preparation). 
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Hence X ~ \/2oi/k ~ 1.5 kpc and the star moves ~ 3kpc 
in each epicycle period. We would not expect spiral arms 
with inter-arm spacing A < 3 kpc to have much effect on 
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